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INTRODUCTION 


The  determination  of  stresses  and  deflections  in  concrete  pavements 
with  joints  and/or  cracks  has  been  a  subject  of  major  concern  for  several 
years.  Since  all  of  the  analytical  (closed  form)  solutions  were  based  on 
an  infinitely  large  slab  with  no,  or  at  most  one  discontinuity,  they  could 
not  be  applied  to  analysis  of  jointed  or  cracked  concrete  slabs  of  finite 
dimensions,  or  with  various  load  transfer  systems  used  at  the  joints  and 
cracks. 

With  the  development  of  high-speed  computers  and  the  powerful  finite- 
element  method,  it  is  possible  to  analyze  concrete  pavements  in  a  more 
realistic  manner.  Various  models  have  been  developed  for  analyzing  pave¬ 
ment  systems  using  the  finite-element  modeling  techniques.  However,  little 
was  done  in  modeling  joints  and/or  cracks  with  load  transfer  systems. 

The  finite-element  computer  program  presented  here  is  based  on  the 
classical  theory  of  meidum-thick  plate  on  Winkler  foundation,  and  is 
capable  of  evaluating  the  structural  response  of  the  concrete  pavement 
system  with  joints  and/or  cracks.  The  model,  which  provides  several  options 
can  be  used  for  analysis  of  a  number  of  problems  which  can  be  summarized  as: 

1.  Jointed  concrete  pavements  with  load  transfer  systems  at  the 
joints. 

2.  Jointed  reinforced  concrete  pavements  with  cracks 
having  reinforcement  steel  at  the  cracks. 

3.  Continuously  reinforced  concrete  pavements. 

4.  Concrete  shoulders  with  or  without  tie  bars. 

5.  Concrete  pavements  with  a  stabilized  base  or  an  overlay,  by 

assuming  either  a  perfect  bond  or  no  bond  between  two  layers. 

6.  Concrete  slabs  of  varying  thicknesses  and  modulus  of  elastici¬ 

ties,  and  subgrades  with  varying  modulus  of  supports. 


I 

GENERAL  ASSUMPTIONS 

The  assumptions  regarding  the  concrete  slab,  stabilized  base,  overlay, 
subgrade,  dowel  bar,  keyway,  and  aggregate  interlock  can  be  briefly  summarized 
as  follows: 

1.  The  smal 1 -deformation  theory  of  an  elastic,  homogenous  medium- 
thick  plate  can  be  employed  for  the  concrete  slab,  stabilized  base 
and  overlay.  Such  a  plate  is  thick  enough  to  carry  transverse 
load  by  flexure,  rather  than  in  plane  force  (as  would  be  the  case 
for  a  thin  membrane),  and  yet  is  not  so  thick  that  transverse  shear 
deformation  becomes  important.  In  this  theory  it  is  assumed  that 
lines  normal  to  the  middle  surface  in  the  undeformed  plate  remain 
straight,  unstretched  and  normal  to  the  middle  surface  in  the 
deformed  plate,  each  lamina  parallel  to  middle  surface  is  in  a 
state  of  plane  stress,  and  no  axial  or  in-plane  shear  stress 
develops  due  to  loading. 

2.  The  subgrade  behaves  as  a  Winkler  foundation. 

3.  In  case  of  a  bonded  stabilized  base  or  overlay,  full  strain 
compatibility  exists  at  the  interface,  or  for  the  unbonded  case 
shear  stresses  at  the  interface  are  neglected. 

4.  Dowel  bars  at  joints  behave  linearly  elastic,  and  are  located 
at  the  neutral  axis  of  the  slab. 

5.  When  aggregate  interlock  or  keyway  is  used  for  load  transfer 
system,  load  is  transferred  from  one  slab  to  an  adjacent  slab  by 
means  of  shear.  However,  with  dowel  bars  some  moment  as  well  as 
shear  may  be  transferred  across  the  joints. 
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MODELING  TECHNIQUE 

For  modeling  the  concrete  pavement  slab,  the  rectangular  plate  element 
with  12  degrees  of  freedom  was  used.  Figure  1-a  shows  that  at  each  node 
there  are  three  displacement  components:  a  vertical  deflection  (W)  in  the 
Z-direction,  a  rotation  (0^)  about  the  X-axis,  and  a  rotation  (8^)  about 
the  Y-axis.  Corresponding  to  these  displacement  components  there  are  three 
force  components:  a  vertical  force  (Pw),  a  couple  about  the  X-axis  ( P0X ) 
and  a  couple  about  the  Y-axis  ( P ©y ) *  respectively.  For  each  element,  these 
forces  and  displacement  can  be  related  by  matrix  notation: 

{P)e=  [Ktop  +  Worn  +  We  ,Dle  <'> 

where  [Ktop]e.  [Kb0ttomV  and  ^Ksub^e  are  the  stiffness  matrices  of  the 
top  layer,  bottom  layer,  and  subgrade,  respectively.  {P}p  is  the  force 
vector  and  {D  ^  is  the  displacement  vector  of  the  slab  element.  For  the 
case  where  two  layers  (slab  and  stabilized  base  or  slab  and  overlay)  are 
bonded,  an  equivalent  layer  based  on  the  transformed  section  concept  is 
used  to  determine  the  location  of  the  neutral  axis  of  the  element.  The 

following  equations  give  the  location  of  neutral  axis  for  bonded  two 
layer  system  using  the  first  moment  of  the  equil valent  area  of  the 


transformed  cross  section. 


2  (ht  +  Vht 


h  .  +  jr-  h. 

t  Et  b 


6  =  I  (ht  +  V  ‘  a  1 

where  a  is  the  distance  from  the  middle  surface  of  the  bottom  layer  to 
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the  neutral  axis,  0  is  the  distance  from  the  middle  surface  of  the  top 
layer  to  the  neutral  axis,  and  ht>  h^,  Et>  and  Eb  are  respectively, 
thickness  of  the  top  layer,  thickness  of  the  bottom  layer,  modulus  of 
elasticity  of  the  top  layer  and  modulus  of  elasticity  of  the  bottom 
layer. 

The  bar  element  as  shown  in  Figure  1-b  with  two  degrees  of  freedom 
per  node  is  used  to  model  dowel  bars  at  joints.  Two  displacement  compo¬ 
nents  at  each  node  are:  a  vertical  displacement  (W)  in  the  Z-direction, 
and  a  rotation  (0y)  about  the  Y-axis.  Corresponding  to  these  two  dis¬ 
placement  components  are  two  forces:  a  vertical  force  (P^)  and  a  couple 
about  the  Y-axis  (P„  ).  The  force-displacement  relation  (including  shear 

8y 

deformation)  for  each  bar  element  can  be  written  as: 


IPlb  ■  ^dowel lb  (D) 


(4) 


where  [Kdowe] lb  is  the  stiffness  matrix  of  the  dowel  bars,  {P>b  is  the 
force  vector,  and  {D)b  is  the  displacement  vector  of  the  bar  element. 

The  relative  deformation  of  the  dowel  bar  and  surrounding  concrete 
is  represented  as  the  stiffness  of  a  vertical  spring  element  (Figure  1-c) 
between  dowel  bar  and  surrounding  concrete  at  the  joint  face.  The  force- 
displacement  relation  for  this  spring  element  can  be  written  as: 


F  =  kdci  a 


(5) 


where  is  the  stiffness  of  the  spring  element  representing  the  dowel - 
concrete  interaction,  F  is  the  shear  force  on  the  dowel  bar,  and  A  is  the 
relative  deformation  of  the  dowel  bar  with  respect  to  the  surrounding 
concrete  at  the  joint  face. 

Neglecting  the  moment  transfer  (if  any)  across  a  joint  or  crack 


wnere  load  transfer  is  achieved  only  by  means  of  aggregate  interlock  or 
keyway,  the  spring  element  show.i  in  Figure  1 -c  with  one  degree  of  freedom 
per  node  is  employed.  The  displacement  component  at  each  node  is  a  verti¬ 
cal  displacement  (W)  in  the  Z-direction,  and  the  corresponding  force 
component  is  a  vertical  force  (P^).  The  force-displacement  relation  for 
a  spring  element  can  be  written  as: 

tP)s  ■  r%9]s <6> 

where  t^gg^s  is  the stiffness  matrix  of  the  spring  element,  and  {P} s  is 
the  force  vector,  and  { D} s  is  the  displacement  vector  for  the  spring  element. 

The  overall  structural  stiffness  matrix  [K]  is  formulated  by  super¬ 
imposing  the  effect  of  individual  element  stiffnesses  using  the  topological 
or  the  element  connecting  properties  of  the  pavement  system.  The  overall 
stiffness  matrix  is  used  to  solve  the  set  of  simultaneous  equations  having 
the  form: 

{P}  =  [K]  {D}  (7) 

where  {P}  is equi valent  nodal  forces  for  a  uniformly  distributed  load  over 
a  rectangular  section  of  the  concrete  slab,  and  { D)  is  the  resultant  nodal 
displacements  for  the  whole  system.  The  generalized  stresses  are  then 


r 
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INPUT  GUIDE  FOR  "ILLI-SLAB"  FINITE-ELEMENT  PROGRAM 

The  "ILLI-SLAB"  finite-element  program  provides  solution  for  the 
deflections  and  stresses  due  to  loading  of  concrete  pavements  vpth  joints 
and/or  cracks.  Longitudinal  and  transverse  joints  may  have  any  or  a 
combination  of  load  transfer  system  such  as  dowel  bars,  aggregate  interlock, 
and  keyways.  ILLI-SLAB  program  is  also  capable  of  handling  a  stabilized 
base  or  an  overlay,  by  assuming  either  a  perfect  bond  or  no  bond  between 
two  layers.  Thickness  of  the  slab,  concrete  modulus  of  elasticity 

and  modulus  of  subgrade  reaction  can  be  varied  from  node  to  node. 

The  wheel  loads  may  be  applied  to  any  slabs,  and  the  stresses  and 

deflections  at  all  of  the  nodes  in  the  slab,  stresses  in  the  stabilized 
base  or  overlay,  vertical  stresses  on  the  subgrade,  and  transferred  load 
by  dowel  bars  are  computed. 

The  concrete  pavement  can  consist  of  1 ,  2,  3,  4,  or  6  slabs 
separated  by  one  longitudinal  and  two  transverse  joints.  The  slabs  are 
numbered  from  1  to  6,  beginning  from  left  to  right  in  the  direction  of 
X-axis,  and  from  bottom  to  top  in  the  direction  of  Y-axis.  Each  slab  is 
divided  into  rectangular  elements  of  various  sizes.  The  elements  and 
nodes  are  numbered  consecutively  from  bottom  to  top  along  the  Y-axis, 
and  from  left  to  right  along  the  X-axis.  Joints  are  treated  as  rectangular 
elements  having  zero  width. 

The  program  can  accept  only  fixed  form  type  of  format  which  is  detailed  below. 
Card  #1 


^  NIX 

N2X 

N3X 

N1Y 

N2Y 

NFOR 

15 

15 

15 

15 

15 

15 
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NIX  =  number  of  nodes  in  X-direction  in  slabs  1  and  4. 

N2X  =  number  of  nodes  in  X-direction  in  slabs  2  and  5. 

N3X  =  number  of  nodes  in  X-direction  in  slabs  3  and  6. 

N1Y  =  number  of  nodes  in  Y-direction  in  slabs  1,2,  and  3. 

N2Y  =  number  of  nodes  in  Y-direction  in  slabs  4,  5,  and  6. 

NFOR  =  number  of  loaded  elements. 

Card  #2  (Use  as  many  as  needed) 

^  XC”(I),~  1  =  1  ,  N1X+N2X+N3X 

8F  10.3 

XC ( I )  =  X  -  Coordinate  of  node  I. 

Card  #3  (Use  as  many  as  needed) 


YC(I) 

Card 


YC( I) ,  1=1,  N1Y+N2Y 
8F  10.3 

=  Y  -  Coordinate  of  node  I. 

#4 


^ NSLAB 

NLAYER 

COMP 

CK 

15 

15 

15 

F10.3 

NSLAB  =  number  of  slabs,  1,  2,  3,  4,  or  6. 
NLAYER  =  number  of  layers,  1,  or  2. 


COMP  =  composite  action  factor,  set  COMP  to  0  if  no  bond 

exists  between  the  slab  and  stabilized  base  or  overlay, 
and  set  COMP  to  1  if  complete  bond. 
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CK  =  set  CK  to  the  value  of  subgrade  modulus  if  subgrade  modulus 


at  all  points  are  equal,  and  set  CK  to  0.0  if  not. 


Card  #5 


BSB 

CEI 

v(D 

F10.3 

E10.3 

F10.3 

CT1  =  set  CT1  to  the  value  of  the  top  layer  thickness  if  thick¬ 
ness  of  the  top  layer  at  all  nodes  are  equal,  and  set  CT1 
to  0.0  if  not. 


CE1  =  set  CE1  to  the  value  of  the  modulus  of  elasticity  of  the 
top  layer  if  it  is  equal  at  all  nodes,  and  set  CE1  to 
0.0  if  not. 

V(  1)  =  Poisson's  Ratio  of  the  top  layer. 

Card  #6  (Read  only  if  CT1  =  0.0,  use  as  many  as  needed) 


T1 ( I ) ,  I  =  1,  ((N1X+N2X+N3X)*  (N1Y+N2Y) ) 
8F10.3 


T1(I)  =  thickness  of  the  top  layer  at  node  I. 
Card  #7  (Read  only  if  CEI  =  0.0,  use  as  many  as  needed) 

^El(I),  I  =  1,  ((N1X+N2X+N3X)*  (N1Y+N2Y)) 

8E10.3 


El (I)  =  modulus  of  elasticity  of  the  top  layer  at  node  I. 
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Card  #8  (Read  only  if  NLAYER  =  2) 


Card  #12  (Read  only  if  N2X  or  N3X  not  equal  to  0) 


^LTDX 

15 

LTDX  =  type  of  load  transfer  in  X-direction,  set 
LTDX  =  0,  if  aggregate  interlock  or  keyway 
LTDX  =  1,  if  dowel  bars 


LTDX  =  2,  if  a  combination  of  dowel  bars  and  aggregate 
interlock  or  dowel  bars  and  keyway. 

Card  #13  (Read  only  if  LTDX  =1  or  2) 


gras 

DOUT 

DE 

DS 

DL 

DJW 

DPR 

FI  0.3 

F10.3 

E10.3 

F10.3 

F10. 3 

FI  0.3 

F10. 3 

DIN  =  inside  diameter  of  the  dowel  bars,  set  DIN  to  0.0  for  round 


bars. 

DOUT  =  outside  diameter  of  the  dowel  bars 
DE  =  modulus  of  elasticity  of  the  dowel  bars 
DS  =  spacing  of  the  dowel ibars 
DL  =  length  of  the  dowel  bars 
DJW  =  joint  width  opening 
DPR  =  Poisson's  ratio  of  the  dowel  bars 
DCI  =  dowel -concrete  interaction 

DCI  for  a  round  steel  dowel  bar  may  be  determined  from  either  Fri berg's  dowel 
analysis  or  from  the  relation  developed  based  upon  a  three-dimensional  dowel 
analysi s: 


a)  Fri berg's  Analysis 


.0.75  n2.5 

DCI  =  - 5-75 - — 

0.041  DU‘ /S  +  0.0004  r-“  w 

b)  Three-Dimensional  Analysis 

E0. 75 

DCI  =  (0.057  -  0.010  D) (0.810  "+  0.01'3T)(T  +  TT/41'4  Wj 

where 

£  =  concrete  modulus  of  elasticity,  psi 
D  =  dowel  diameter,  in. 
h  =  slab  thickness,  in. 

W  =  joint  width  opening,  in. 

K  =  modulus  of  dowel  support,  pci 

Card  #14  (Read  only  if  LTDX  =  0  or  2) 


HU 

E10.3 

Agg.X  =  aggregate  interlock  factor.  Use  a  large  value,  i.e., 

g 

Agg.X  =  10  for  keyways. 


Card  #15  (Read  only  if  N2Y  not  equal  to  0) 


03 

15 

LTDY  =  type  of  load  transfer  system  in  Y-direction,  set 
LT0Y  =  0,  if  aggregate  interlock  or  keyway 
LTDY  =  1 ,  i f  dowel  bars 

LTDY  =  2, if  a  combination  of  dowel  bars  and  aggregate 
interlock  or  dowel  bars  and  keyway 
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Agg,Y  =  aggregate  interlock  factor.  Use  a  large  value,  i.e., 

O 

Agg,Y  =  10  for  keyway. 


Card  #18  (Read  NFOR  times) 


mm 

PRS 

XI 

X2 

Y1 

Y2 

■ 

15 

F10. 3 

F10.3 

F10.3 

F10.3 

F10.3 

■ 

NEL  =  element  number  of  the  loaded  element 
PRS  =  tire  pressure 


XI ,  X2  =  lower  and  upper  limits  of  the  loaded  area  in  X-di recti  on, 
in  element  coordinate  system  (local  coordinates). 

Y1 ,  Y2  =  lower  and  upper  limits  of  the  loaded  area  in  Y-direction, 
in  element  coordinate  system  (local  coordinates). 
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TYPICAL  EXAMPLES 

Five  example  problems  are  presented  here  to  illustrate  various 
applications  of  the  "1LLI-SLAB"  finite-element  program,  for  determining 
the  response  of  the  concrete  pavement  system  due  to  loading.  Following 
each  example  problem  the  finite-element  mesh  used  for  each  problem,  computer 
inputs  and  outputs  and  a  summary  of  maximum  stresses  and  deflections  of  the 
slabs,  as  well  as  load  transferring  systems  are  given. 


Example  1,  Jointed  Concrete  Pavement  with  Dowel  Bars 


Slabs  -  Two  25  ft.  square  panels. 

Thickness  =  12  in. 

Modulus  of  elasticity  =  5  x  10®  psi 
Poisson's  ratio  =  0.15 

Load  Transferring  System  =  Round  steel  dowel  bars 
Inside  diameter  =  0.0 
Outside  diameter  =  1  1/4  in. 

Spacing  =  15  in. 

Length  =  24  in. 

Modulus  of  elasticity  =  29  x  10®  psi 
Poisson's  ratio  =0.29 
Dowel -concrete  interaction  *  2.4  x  10® 
Joint  Width  *  0.10  in. 

Modulus  of  Subgrade  Reaction  =  200  psi 
Load  =  50,000  #  at  the  center  of  joint  in  Slab  1 
Contact  pressure  =  222  psi 
Contact  area  *  15  ft.  square 
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- 0.-44-024-4*4-- 

0. 840269E 

-02- 

02 

12. COO 

c.o 

0.  0 

12. 000 

c.  c 

0.  0 

O.C 

c.o 

0.  47  50152 

0  3 

12. 000 

0.0 

-  0.  475015E 

03 

12.000 


-0--3-55-80UE — 02- 
■0. 355808E  02 
0.  0 
0.0 

0. 502431E-03 
0. 50243  IE-03 
0.  0 


0.0 

12.000 

..Jt2^CUUL 

12.000 

0.0 

-1-2 . 00  0  - 
12.000 
12.000 


■  0. 84C205E  02 
0. 64  020  52  02 

-0-0 - 

0.0 

■0.501953E  02 


-0. 355812E  02 
0.355B12E  02 

0.0 _ 

0.  0 

■0. 293333E  02 


12.000 

12.000 

12.000 

0.0 
12. 000 
12.000 


0.  C 

12.000 

12.-000 

12.000 

-12-J-OOO- 
12.000 
12. COO 


12.000 
12.000 
12. 000 

0.0 
12. 000 
12.000 


-0. 622276E  02 
0. 622276E  02 
- 

0.0 

•0.  217354S  01 


-0. 2B6678E  03 
0.  0 
0.  0 

•0.217074E  01 
0.217074E  01 
0.0 


-0. 472234E  02 
0.472234E  02 


-0-i-3-6-54-84-E— 03- 
0. 365434E-03 
0.0 
0.  0 

0. 842369E  02 
-0. 842369E  02 
0.  0 


0.0 

12-000 


■0.622232E  02 
0. 622232S  02 


0. 472233E  02 
-0 . 472233E  02 


12.000 
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12. 000 
1? . 000 
12. COO 


0.  0 
12.000 
12.000 


0.  1913581  02 
0.0 
0.0 


-0.317242E  02 
0.317242E  02 
C.O 


■0. 154169E  02 
0.  1 5  4  1 6  9  E  02 
0.  0 


0.  193660E  02 
0.  0 
0.  0 


-0.362705E  02 
0. 362705E  02 
0.0 


o.c 

-0 . 6894  6  1  E 

02 

0.  163552S 

02 

-  0. 406333E 

02 

12.000 

0.689461L 

02 

-0. 163552E 

02 

0. 486333E 

02 

12.000 
0.  C 


12. 000 
12.000 
12. OCO 


0.0 

12.000 

12.000 


0.0 

•  C.  108429  E  03 


0.689456E  02 
0.0 
0.  0 


•C. 3172 50 E  02 
0.317250E  02 
C.O 


0.  0 

0.  525973 E  02 


-0.  163574E  02 
0.  0 
0.  0 


■0. 154159E  02 
0.  154159E  02 
0.  0 


0.0 

•9.237  69  5E- 


-0. 486337E  02 
0.  0 
0.0 


0.362710E  02 
-0. 36271 OE  02 
0.0 


0.0 
i2. ooc 
12.000 
12. COG 


-0.  1 9 1 367 E  02 
C.  191 367  E  02 

do 

-Q.128979E  02 


0.198666E  02 
-0. 198666E  02 
0.0 
0.0 

-0. 145887F  02 


12. 000 
12.000 
12. COO 


0.0 
12. 000 
12. 000 


0 . 288563E  02 

0.0 

0.0 


■0.511912E  02 
0.511912E  02 
0.  0 


■0.6881B6E  00 
0.0 
0.  0 


0.  499480E  01 
•0. 499480E  01 
0.0 


0.1 92242E  02 
0.0 
0.0 


-0.  105093F  02 
0.1 85093E  02 
0.0 


0.0 

1.?.  coo 

H2-r000- 

12.000 


12.000 

12. COO 


•  0. 594073  E  02 
0.5948731  02 

-G-r0 - 

C.  0 

-C. 5119 23 E  02 


0.8253752  01 
■0.  8  2  537  51  01 

-Or-e - 

o.  o 

0.  49  9  48  3  E  01 


0.1051092-03 
■0.  1 05 1 09  E- 0  3 
- 

0.0 

0. 1 B5094E  02 


12. OOC 

0.2 88568 E  02 

-0.6885432  00 

-0 

. 1922442  02 

12. 000 

0.0 

0.  0 

0 

.0 

12.000 

0.0 

0.0 

0 

.  0 

0.  14588BE 

02 

-0. 145888E 
0.0 

- ft-,  A . 

02 

7  1 

0.  0 

-0.241 105E 

01 

0.  0 

-0.919507E  01 

-  ! 

wsmmmm 

0.^411051 

01 

0.0 
o  5 

0.919507E  01 

Q  Q _ 

ft 

12.000 

c.'o 

o’.  6 

6.' 6 

72 

0.0 
i  o  nnn 

-0.  1  195031 

02 
r\  a 

0. 4? 5309E  01 

r\  n  -i 

- 0. 795264 E  01 

A  ni 

J  —  -  - 

12.000 

c.o 

0.0 

0.0 

_ 

12.000 

or 

0.0 

-r 

02 

0.  0 

n  iiq^qt  _nj  _ 

0.0 

-O  L  ^ 17  22  3E  01 

12l 000 

0. 203039E 

02 

-0.  1 1  9208E  01 

0.51 7223E  01 

12. CC0 

0.0 

0.0 

0.0 

12.000 

0.0 

0.  0 

0.0 

74 

0.0 

-0.236805E 

02 

0.8141122-02 

-0. 139564E-04 

12. COO 

0.236805E 

02 

-0. 8141122-02 

0. 139564E-04 

12.000 
ii  r  nn 

C.O 

0  0 

0.0 

Q  -Q. 

0.  0 

0  u  o  - 

75 

0.  0 

-0.  203039E 

02 

0.  119224E  01 

0. 517216 E  01 

12.000 
...  -  -p  rn/v 

0.203039E 
0  0 

02 

-0.  1  19224E  01 

0  0 

-0.517216E  01 

0  _  Q _ 

12.000 

0.0 

6."  6 

6.5 

1  76 

0.  C 

-C. 1 19505E 

02 

0.435313E  01 

0. 795269E  01 

- 3  "3 - 1*  v  - 

12. COO 

- WTJWmTg- 

-0.306549E 

- 

01 

- Vm  -v - 

0.  0 

- it  V  <7  J  H  n 

0.1 70084E 

01 

12.000 

0.0 

0.0 

0.0 

12.000 

0.0 

0.0 

0.0 

B6  0.0 

0.1 39047L 

01 

0.  16  99  0  IE 

01 

—0.71 4609E 

00 

12.000 

-0. 139047E 

01 

-0. 169901E 

01 

0.7  14609E 

00 

12.000 
“t  A  AAA 

0.0 

0.  0 

A  ,  A _ 

0.0 
o  n 

3  7  0.0 

C. 9191370 

00 

0.527423E 

00 

-0.  145617E 

00 

12.000 

- lj—33-0 — 

-0 .  j)1 9 1 37E 

00 

-0. 527423E 

00 

0.  145617E 
- 0^0 - 

00 

- ^ - 

- 0-.-0 - 

12.000 

- 0.919  139-fc- 

-0.9 19 119 r 

-aa— 

00 

— a  ^4  a  an  3i- 
-0. 527403;; 

-90— 

00 

- 0^Ui-5<.a-J7 

-0.  1 45603  E 

-ou - 

00 

12.000 

C.O 

0.0 

0.  0 

12.000 

0.0 

0.0 

0.0 

> 

90 

0.  C 

C.  1  39047  E 

01 

0.  169899E 

01 

0.7  14596E 

00 

12.000 

-  0.  1  39047E 

01 

-0.  169899  E 

01 

-0.714596E 

00 

.  \hm 

8:8. 

8:8 

8:8 

91 

0.0 

C. 306539  E 

01 

0.0 

0.  170082E 

01 

12. 000 
- 12-.- 896 

- C. 306  53  9  F 

- 9rfr - 

01 

0.  0 

- £h-9 - 

-0. 170082E 
- frr-0 - 

01 

92 


12. COO 

0.  0 
12-.-&0C 


0.0 

C.O 

-O.-fr- 


0.  0 

0.  0 
-fr*-fr- 


0.0 

0.0 


-9-3- 


12.000 

12.000 

-0-.fr 


C.O 

0.  0 

-0-5“0- 


0.  0 
0.0 

-0.-689923?  -08- 


0.0 
0.  0 

-9.^9832-95-04- 


12. 000 
12.000 
12.000 


0.0 

C.O 

0.0 


0.  689923E  00 
0.0 
0.  0 


0.258320E  00 
0.  0 
0.0 


94 


0.0 

12.000 

12.000 

-^K'-r-frOfr- 


0.0 

0.0 

0.0 


0.  79442  IE-01 
-0. 79442 1E-0 1 
0.  0 

- 


-0. 597485E-01 
0.597485E-01 
0.0 

_0_0 - 


95 


0.0 

1g. COO 


0.0 

0.0 

-frvfr- 

c.o 

C.O 

-fr^fr- 


0.  26124BE  00 
0.  26  1248E  00 
-frv-e - 


-0.224948E-05 
0. 224948 E- 05 

-frv-9 - 

0.0 

0. 597552E-0 1 
-0.-597552F— 01 


96 


12.000 

0.0 

1 2.-600 


0.  0 

0.  79415 6E-0 1 
-9.704  166  £-04- 


12.000 

12.000 


C.O 

0.0 


0.0 

0.0 


0.0 

0.0 


- - -j-j - 

-  V  A  V 

12.000 

12.000 

12.000 

- - 

0.0 

0.0 

0.0 

- - i  r — irv — 

-0.689901E  00 

0.  0 

0.  0 

- vne-vO-a-i-ot — - 

-0.258318E  00 

0.0 

0.0 

98 

0.0 

C.  0 

0.0 

0.0 

12.000 

0.0 

0.  0 

0.0 

12-000 

C.O 

0.0 

0. 0 

- 12*930 — 

— - 

- 0^-9 - 

— o^o - 
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NO  Dr 


TR  A  KSFI’FRL  E  LOAD 


NODE 


TRANSFERRED  LOAD 


Example  2,  Jointed  Concrete  Pavement  with  Aggregate-Interlock 
Slab  =  Two  25  ft.  square  panels. 

Thickness  =  12  in. 

Modulus  of  elasticity  =  5  x  10®  psi 
Poisson's  Ratio  =0.15 

Load  Transferring  System  =  aggregate-interlock 

Aggregate  Interlock  Factor  (Agg)  =  5  x  10^  psi 
Modulus  of  Subgrade  Reaction  =  200  pci 
Load  =  50,000  #  at  the  center  of  joint  in  Slab  1 


Contact  pressure  =  222  psi 
Contact  area  =  15  in.  square 
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STRESS 


Example  3,  Jointed  Concrete  Pavement  with  Keyway 
Slabs  =  Two  25  ft.  square  panels 
Thickness  =  12  in. 

Modulus  of  elasticity  =  5  x  10^  psi 
Poisson's  Ratio  =  0.15 
Load  Transferring  System  =  keyway 
Modulus  of  Subgrade  Reaction  =  200  pci 

Load  =  50,000  #  at  the  center  of  joint  in  Slab  1 
Contact  pressure  =  222  psi 
Contact  area  =  15  in.  square 


37 


G* 28  7327-6  U/M  OK 

FORTRAN  Coding  Form  Printed  in  U  S. A 


RESS 


Figure  8.  DISTRIBUTION  OF  STRESSES  ALONG  THE  KEYED  JOINT 


Example  4,  Keel -Section  for  Runway 
Slab  -  25  ft.  square  panel 

Thickness  of  the  center  =  12  in. 
Thickness  of  edge  =  18  in. 

Modulus  of  elasticity  =  5  x  10^  psi 
Poisson's  ratio  =  0.15 
Modulus  of  Subgrade  Reaction  =  200  pci 
Load  =  50,000  #  at  the  center  of  thickened  edge 
Contact  pressure  =  222  psi 
Contact  area  =  15  in.  square 
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Figure  10.  FINITE-ELEMENT  MESH  CONFIGURATION  USED  FOR  EXAMPLE 
PROBLEM  4 


Example  5,  Stabilized  Base 


Slab  =  25  ft.  square  panel 
Thickness  =  12  in. 

Modulus  of  elasticity  =  5  x  10^  psi 
Poisson's  ratio  =  0.15 
Base  =  Cement  stabilized 
Thickness  =  6  in. 

Modulus  of  elasticity  =  1  x  10^  psi 
Poisson's  ratio  =  0.25 
Composite  Action  Factor  (COMP)  =  1 
Modulus  of  Subgrade  Reaction  =  200  pci 
Load  =  50,000  #  at  center  of  edge 
Contact  pressure  =  222  psi 
Contact  area  -  15  in.  square 


46 


STRESS 


PROGRAM  LISTING 
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5X7FKNA1  PPOGPr 
Villr  ((.97) 

- ■«,  7-  f o 3P  A  " — (/VirCt; 

1  //.:>cx,** 

2  //, 2  CX , ' * 

3  //,?0X,'* 

-* - //T2«*fM. 


5  //.  2  CX  , 

RI  AD  <5,?-00)  N  1 X,  F:2X  ,  F 3X  ,  K  1Y  ,  N2Y  ,  NFOR 
c.00  FORMAT  (feI5| 


JOINT  PROGRAM 

FINITE  E1EKEKT  ANALYSIS  CF  CCNC.  PAV. 
AMIF  H.  T  AEATA3 AI l 

-iPNl  Vf  SSI-T  Y~  <HF— IrfctiHC  IS — - V97-7- 


*  • 

*  •  ‘ 
*  •  ' 


-*»  x -  f  I  x-»r 


r3"X- 


NY=  K 1 Y  +  N  2Y 
IDM9  =  i:x*FY 
IDN8=IDNS*3 

irM6=icre*iEN7 
IDM5  = (NX-1) *  (NY-1) 

xiau=yycp  _  _  ^ 

iCAtL  ccf.rz^  (incGnn,iM | ^ i E5 6 * f 5 n ” i f) m 4 ? idm^i d a « , 

1  iDr«,iDM9, 

2  I P M 3, 1 DM9* 3. IDM9* 3, I  DM  9* 5. 1 DM9 *5,  ID  19*  5.  I DM9* 

-3 - I  D  f'9  ,  IDil9  ,  IP  t*.9  ,  NX ,  r! V;-9-^H-X .  0  ,  N-S^-yO  ,  H39C  ,-Q7R  1  Yt 

«  0, KFCK,0, NX, 0,IDM9, 0, I DM3 , 0 , 1 E«7 , 0 , 1 D  M  6  ,  0, IDS 5,0, ID. 14) 

STOP 
END 


52 


inp 


2  XC,yC.];iXtl)2X,lljXfKli,H2Y,NY,NFOR,TlX,IDB9,lDH8,IDtl7, 


t. - I  BMP, IE 

INTEGER  CCMP 
DIMrl.’SICN  X  (3DM9)  ,Y  (IDM9) 
1  F2  (ILM9)  ,SUb  (ID 

■? - sr-frew  - 

3  Y  1  |1DM4 


4 imf 


,  V  ( 2)  ,T1  (IEM9) 
.  I'9)  ,  U  (2)  ,R  (2) 

TDItfr'  ‘ 


,T2  (IDM9 
N  (9)  ,  SF  (  .  . 

[-MOi'.XIl 


STFT1  (IDM9 
TT  (IEM9)  ,X~ 
-YihtttMl 


rsv-iiim6) ,  hel  (±pn«rvP«3jr&fl<n  .  x  wi-ginttr — r - 

,FC  jIDM9K  P(IDM8)  ,  STR 1  (IDM9  ,  M,  STR2JIDH9, 3) , 

M9 , 5)  . STRBl JIDM9,  5}  ,STKT2 (IDM9,  5)  . STBB2  (IDH9,5)  , 
F(IDH9)  ,  YF  (IDH9)  ,XC  (NX)  ,  YC  (KY)  ,XZ  (IDM4)  , 


2  FORMAT  (EF  10.3) 

3  FORMAT  (//, 1  OX , ’ NO.  OF 

*  nxfM!0.  CF  NODES 


-r-wr, 

*  10X. 

♦  10X. 
4  FORMAT 


NODES  IN  X-DI R ECTION 
IN  X  DIRECTION  SLABS 


SLABS  1,4=', 15,//, 

hZzUi&Si _ 


//7 


NC.*  C?  NODES  IN  Y"DIRECTIC’!  SLABS  1,2 '3=  *'1*5, // 
NO.  CF  NODES  IK  Y— DIRECTION  SLABS  4, 5, 6=*, 15) 

,  1 CX,  ’ X-COOPBI NATES  ARE: • ,//, 1C (3X, FI 5. 3  ‘ 

'  «  '  a-,.  mm-  ** _ »  r,n  .  *  w  .  .  r  4  r\  1>  w  m  a  — 


5  FORMAT 


1  OX  .  1  Y-COORDI  KATEg-Afr£t 1  ir/H 


Tn.  luxrx 

8  FCFMAT  (315 ,110.3) 

9  FORMAT  (//. 1 CX, *  NO.  OF  SLABS=’ 

*  //. IPX, rCCME.  ACT IC  N= 1  .15) 

4 Ofn)i  >  A 4-A-V »-T>  r  ftTW-g — fVF — 


u- 


12 
13 
16 


°CR<‘f4-t//.  1  OX.  1  P80P9RTE^B 
FORMAT  (F10.3,F10.3,F10.3) 
FORMAT  (/, 10X, ’ P01SSCN  RAT 


£&■ 


,  I5,//,'10X,  *  NO.  CF  LAYERS=  •  ,15  , 
THE  TOP-  LA  Y'ER— i~S:-*  ) - 


FORMAT  (/, 10X, ’THICKNESS 


18  FORM  A* 

19  FORMAT 

*  8(15,110.3)) 

22  FORMAT  (/,1CX,* 


i/  t  > «  *  , 

l/,iox,* 


10  OF  TOP  IAY£R=’,F10. 
CF  TOP  LAYER=* , F10.  3) 


-2  4 -FORT*" 


THICKKISS  CF  TCP  LAYER  AT  THE  NODES  IS:’,/, 
MODULUS  CF  TOP  LAYEE=*  ,E1 0.  3) 


■-  tBEIO. 

25  FORMAT  I/,10X,' EODULUS  OF  TOP  LAYER  AT  THE  NODES  IS:’,/, 

*  8(15, riO. 3)) 

27  FORMAT  (//, 1 CX. ’PROPERTIES  OF  THE  BOTTOM  LAYER  IS:’) 

29— •  (/,  TO?.,  PCIOSCK"  RATIO  -€P  DOTTCM  LAYER  I  St»  ,F-10.3) 

31  FORMAT  -  -  “  - - -  " 

33  FORMAT 

* 

-96 


[/,  10 X,  ’THICKNESS 
I/, 10X, 'THICKNESS 


CF  BOTTCM  LA  YER= ' , FI  0. 3  I 
CF  DOTTCM  LAYER  AT  THE  NODES  IS:’,/ 


8(15,110.3)1 

OP  MAT  (/,  1 OX ,  MODULI' 3 — OF  -BOTTOM -LA-YE R-  1  ,E1-0.  3 
CFMAT  (/, 10X, 1  MODULUS  OF  BOTTCM  LAYER  AT  THE 


38  FCFMAT  </, 1  OX, 1  MODULUS  CF 
*  8(35,510.3)) 

42  FORMAT  (/,10X,‘ SUEGRAIE 
MM-FO-FtM-AT — |8FK%-3| - 


BOTTCM  LAYER  Ai^THE^ODE S  IS:’ 
BCDU1US=' ,F10.3) 


/. 


45  FORMAT  (//,1CX, 'SUBGRADS  MODULUS  AT  THE  NODES  IS:’,/, 
8  (I5.F10. 31) 

ID AT  (2F10.3,E10.3,4F10.3,E10.3) 


6  FORK 


■9  -PORH-AT  (//>1 OX? ’PROPERTIES  OF— THE  DO«ti  EARS- ABB;  * 


•INSIDF  D1A.  =  ’.F10.  3./,  10X,  'OUTSIDE  DIA.  =  •  .  F  1  0.  5, /,  1  OX  ,  „ 

•  rCEULUS  OF  ELASTICITY’, £10. 3, /.  10X,  ’SPACING=’ ,FT0. 3,/, 10X. 
’L1NGTB=’ , F 10. 3, /X1 OX, ’POISSON  R AT  10= • XF 1 0. 3, 


T  Ll^win-  ,  .r  iv.jf/,  <va,’POISSOH  i1.  n  *  .  w  -  ,  *  i 

* - /,10V,  »  DOWEL— CCRC  PETE  INTE-ft*C'HOW=-t  x  &-10.-3 


*  //,1CX 
6  1  FOP  ■- 
Ii2  FOR 


' J Cl NT  KIDTH=’ , F 10. 3) 


/  /  |  I  VAf  U  VA  I*  HAL/ 

MAT  (15,5110.3) 
MAT  (//.  13X,  ’  ELEI1 
-  li,  i.V4-rnn 


I/  /  I  I  J  A  I  (  -A 

21 1  »  ¥  1  -  C  OO  R-  t-J~r‘  X>- 
(/,10Xf15,5X.5F10.  3) 
/>,1 1*,  ’flOD^.SX,  >b 


ENT'  5X,  ’ PRESS'  ,2X, ’X1-C00B. • ,2X, ’X2-C00B. 

P-COOR*  «-) - — - 


64  FORMAT  1/  ,  ■  w  a  ,  )  j  a.  j;  iw.  jj 
70  FORMAT  (//, 1  1  fc,' NODS f , 5X ,’ DEFLECTION* , 5X ,' X  ROTATION’, 
*  5>:,’y_ROTATlCN’ ,  1  5X, '  SUBGRADE  STRESS',/)"  _ 
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wri. 

GC  TO  37 

-y.— E^-a-9-h»  1  ,  ID  K  5 - 

39  T2  (I )  ■=?. 

37  IF  (CK.1C.0.)  GO  TO  40 
DO  41  I=1,IDE9 
~SffB -fTJ 


trt- 


rTi-^rn- 


WRITE  (6,42)  CK 
GC  TO  43 

40  RFAD  (5,44^  (S OB  (I )  ,  1=  1 ,  I ET9 ) 

- WftTf-' (t,45)  -(iTjtfSil) ^1  =  1ri-£^9)- 

43  CONTINUE 

IF  (NS1AB.EQ.  1)  GC  TC  49 
IF  ^’2;.rg.°)  GO  TO  999 


prat— f5r 


5t?o01  IT 

MK 
} 


tnr 


WPITE  J  . 

IF  (LTCX.E0.1)  WPITE  (6,2  12 
IF  (ITEX.EQ.Q)  WRITE  (6, 4 11 


tr-r* 


WRIT’ 


hrrrT^ 


IF  (1TDX.EQ.0)  GO  TO  201 

P.JAD  (5,6)  DIV'X,D0UTX,DE.X,DSX,D1X,DJWX,DPRX,DCIX 

I  '  ' 


_  I  NX^POUTX,  DEXjlDSX,DLX,EPRX,ECIX,EJWX 
^  +  ft-efiTX*  +  4  DTHf»*S) - - - - - - - 


IF  (EINX.iC.  -)  GO  TO  149 
SHAX=C. 35061 1*  (DOUTX-DIKX) * ( (DOUTX**2-fDINX**2) **2) / ( (DOUTX* 
1  El  ICX)  ♦♦  3) 

QO  TC  150 


SHAX=C.70685£*  (DOUTX**2^ 


149  .  _  _ 

150  ?  YX=  24. *EIX*  (1.+DPRX)/  (SHAX*EJWX**2) 
IF  (LTEX.FQ. 1)  GO  TO  999 

— fiV  9ftZ> — S- 


-20T  PFftt  ~t5;20 


WRITE  (1 , 203)  AGGX 
999  IF  (K2Y.FG.0)  C-0  TO  49 
REAE  (5.200)  ITDY 
-  & 


-WPTT 
IF  ( 
IF 
IF 
"K 


(ft. 191) 
LTEY.FQ.  1 
LTEY.EQ. C 
2 
-6 


WRITE 

6.2121 

i  W  FIT’S 

WPITE 

- fifS  *»D - 

mm 

LTDY.FQ. 
lx.-TCY.EQ. 


READ  (5.6)  DlJiY.DOUTY.Di Y . DS Y, DLY. DJWY. DPR! . DCIY 
WRITE  It, 7)  EIKY,DOUTY,DEY,D5Y, DLY, DPBY, DCIY, DJWY 
EJY=0.C4S087* (EOUTY**4-DINY**4) 

— rr-fftTTYT-Eft:  eo — po  tc  349 — — - 

SHAY=  C. 3506 1  7* (DOUTY-EINY)* { (DOUTY** 2 *DI NY**2) **2) / {  (DOOT  + 
*  El  RY) ** 3) 

rr  ic.n 


GO  TC  350 

349  SHA¥-€a70685€* t9WTY»»2) - 

350  EYY=:24.  *EI  Y*  (1.  +DPBY) /  {SHAY*DJWY**2) 


55 


iP  rnpy, F£.  n  go  to  49 
301  nr/.r  (5,;o2i  Kggy 

- -WF3TE — (Ft2-OS)— A-GGY - 

44  U  (1)  =0.r<  (1.-V  (1)  ) 

II  (2)  =  C.C*(1.-V  2)  f 
DC  to  3= 1.2 


54  YF(I)  =  C. 

DO  55  1=1. IDEA 

55  P  (I)=0. 

- C~ft£r~£'Ti  ? F'  t~3  D fr5 .  H  Xt  N' 3  Y  7*’ 2  Y7-I6Y-7 - 

*  X.Y.IMO ,NSLAB,N1X,H2X,N3X,SF,SUB,T1,E1, V, 

1  T2.£2,U,R,S,N,IDM7.ST,IDM6, COME, NL A YER , 

2  EIJIX.rOLTX.DEX,  DSX.DLX,  LJWX.DIX,  DPRX.FYX.LTDX,  AGGX  . 

- 3 ~TT~?4  Y  rD  O  t T  YtD  rY  ,-D  ST  rB  L  Y  rE  J  w  ytbi  Yt  D  P  R  Y  v?  Y  Y  7E  T  C  Y  V  *  G  G  Y^  Cl  x^-d  Cl  Y )  - 

CALL  UPTEI  (ID!i7,IDH8,S0,£l,IDM6J 
DC  6  C  1= 1 ,  I E  F  4 

60  READ  (5.611  » 1 1  (I )  ,  PBS  (I )  ,  X  1  (I )  ,  X2  (I)  ,  Y  1  (I )  ,  Y2  (I)  * 


DO  63  1=  1 ,  ID  F  4 

63  WRITE  (6,64)  }  £ L  ( I )  .  PES  (I )  ,  X  1  (I)  ,  X  2  (I )  ,  Y  1  (I )  ,  Y  2  (I ) 

CALL  1CAE  (IEPI4.NEL.  FRS ,X  1, Y 1, KY ,X, Y , IDE9,  F.IDB 

- C* L T.  rlTTP  tT-t-B^TT  DaVl  .T  .  SD^^TrlDHtr) - — 1 - 

WRIT’-  J6,76l 
DO  90  1=  1,  ID*9 
FC  (I1=P  (  (1-1  )  *3*U*S3B  (I) 


-JJ  S'  C.  d.  f  A  .XU.1U.A#!  f  •  ,  J,  '  f 

»:si  AB.  X  .  YF, DEX ,ESX  ,  CIX , DJ WX , DPR X.F YX, LTEX , AGGX , 
rEY.ESY.flY, D  JV  Y  , DPRY.FYY ,  LTDY , A  GG  Y , PCI X , DCI Y) 


DO  80  . 
WRIT  I 
WRITE 


(STRT1  (I.J)  ,J=1.5) 

.  (STFB1  (I.J)  .J-1,5) 


WRITE  (6,190) 
WRITE  (6,83) 

J1  1=  (SIX-1)  *FY*1 


IF  (N2Y.IQ.C1  GO  TO  9C 
IF  (ITEY.  ?Q.  t)  GO  TO  90 
WRITE  (  G> ,  1 9  1 ) 

- ttfrlTE-  |f.G3) - 

nx=kix+k2x+n3x 

JJ1  1«  (HX-1) *  Y«N1Y 

00  8fif  K=»!1  Y  ,  JJ  1 1  ,  NY 

m=K 

LL2=K+1 

P8C  W?  IT  I  (f ,  85)  Lll.YF(lI.I)  ,LL2,YF(LL2) 
90  CONTINUE 

- irtTory - 

EKE 
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r.UDKCUTJKS  SHF?  (IDG5,  NX  .HI  Y,  K2Y,  IIY. 

„  X, f,IDR$,R§LAB,NlX, N2X.N3X.SF, sup, 

.  T1, 21^72,1  2, ll.r:.S,K,IDM7.  ST, IPn6,C0nP,NLAY£R, 

- Pi  NX,  TrP-  WT-frfPf-X-7-P  r.  X  -,-PL  X-i-C  J+i  X-r&5-X-rBPr-XTi!i  X  ,3.-5  D  X  ,  A  GGX-, - - - 

CSY.D1Y,  C.7WY,DI  Y  ,  DP?  Y,  FYY  ,  LTt  Y ,  AGGY,  DCIX.DCIY) 

lHTtuzi*  iCnP 

DIHFNSICJ!  X  (3tr9)i  ,Y  (IEK91,SF  (90)  ,SUD(IDK9)  ,T1  (IDK9)  ,E1  (IDH9)  , 

Dr~10~I  =  1  j  ,n  (-2M  (jo,2)  ^-f4)-rST44-6«H— 

NT  1=1*  (I -if/  (NY-1) 

NI 2=  NI  i  +  1 


X  (NE3 


IF  (B.FG.O.)  .  GC  TO  601 
GO  TO  1C2 

IT  (HEX. iC*  0)  GO  TO  510 


1)  DJWX=0. 


GO  TO  101 


IF  (DJHX.LF.C. 

DCX=DCIX*B/D  X 

Et’X=  1 2.  +  E  PX/ (D JWX**3  *  (1.  +FYX)  ) 

C2x=6.*rrx/ (tjvx**2*  (I.+fyx) ) 

E3X=  (4.+FYX)  *rEX/(DJWX*(1.+FYX 
D4X=  (2.-FYX) *EEX/ (DJWX* ] 1.+FYX 
C1X=  (1./  (  (1./ECX)  +  (1  ./E5X))  ) 


SF  (3) =E2X 
SF 
SF 


=  -E2X 
=  E4X 
=  E1X 


=  E4X 
=  E1X 
=-E2X 


GGX*  B 

-ST  (  1)  ♦  AuX 
)=S F ( 19  J-AGX 
)=3F (37) *A OX 
SF  (55) =SF (55J-AGX 
S5  (6i|)=SF  (64)  «AGX 
SF  (6 2)  =  S F  (82)  *AGX 


S  (21  ,0 
S (22 ,0 
S  (23 , 0 
S  (2«,0 
S 

S  (27,0 


♦  (13.*BS-8.*U  (J)  )  *A S** 

♦  (— 60.*BS+30.*AS-30. *V  (0)  -8**.*U  (0)  ) 


*  (30.  +BS  +  6.*U  (J)  )  *A2 

1.0} 

*  ( 10.  *AS-8.*U  (0)  )  *ES4 

M,0) 

*  (10.»BS-2.*U  (0)  )  ♦AS'* 


H7TURK 

f»-’D 
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■■pm 


, su.st,2p:i6) 


i 1=1 r*7- i 
2?=irrc-iL^7+ 1 
DC  1000  2=1,32 
-3=-rT^r*irr;7  «-t 
SU  (J )  =sT  (J) 

cc  icci  k=i,;i 

SU  ( J  +■  K)  =  ST  {J  -»K) 


?~>K)-ST-(-J«H1/SU  (.1) 


-toot  i>~i .. 

13=1  1 
DC  1000  1=1.21 
T=J+1*II:  07-1 


1002 

1000 


-&C-  1-0C2  ~K=  1,33 


ST  (M  +  K)  =  ST  (tt  +K)  -ST  (J  +  K  +  2  1-23 }*SU  (J+L) 
23=33-1 
22=12+1 
-3-4~I-l  1 


DO  1003  7=22 , 2E.,18 
J=  {I-  1)  *IEP7  +  1 
SU (J)=ST  (J) 

-frF-fF  14)  13C4,  1004,  1C05 


DC  0  06  K=  1  .  j  *4 
SU  (J  +  K)=ST  (0  +  K) 

CC  1007  K= 1,11 

T  (0  *  K)  *S?-f-3-rKj-/SO  ( Jf 


1004 
1006 

1005 
1  C-9-7 


13=1 1-1 
IF  {I— I £  :18) 
1008  DO  1010  L= 1 


~=hA 


1CC0, 1003, 1006 
34 


1011 


DO  1C  11  K=1, 3  3 

ST  (M  +  Kl^ST  (A  KJ-ST  (J  ♦  K  +  1 1-23-  1)  *  SU  (0  +  L) 
tF  1 1  13) — FfrIO, -10  10,  1012 - 


1010 

1012 

1003 


CCNTINUF 

14=24-1 

COOT  2  1JUF 
- 


EJ’D 
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- :  Hi  IM?*,!  - 

’»e  20C6 

AF1=ME  (1)  +  (Ml  (I) -1 )  /  (KY-1) 

- X  »2="*F  1«-1 - 

Nl3=»F1+NY  , 

KE4=  M  1  4  NY  +  1  * 

_ *f  _ 

^pAsCI 

F2=PRS  (3) 
F  ?  =  P  RS  1) 

— r 

*  ( 
*  ( 

t'iu  |-£V(I))*  (Y2  (I)  -Y  1  (I)) 

X2  (3)  **2-X1  (I)  **2 )  *  (Y2  (I)  -  Y  1  (I)  1  /2. 

JUi) L£L*A2 

*5=FFS  I  M 
F6=P  RS  (I) *  | 
*7=TKS  (3)  *  | 

X2  I 
X2  (I) 

illil 

*  *  2-  X  1  I  *  *  2)  *(  Y2  ( X)  *  *  2-Yl  ( I T*  *  2 )  /«. 

-XI  (3)  )  *  (Y2  (I)  **3-Y1  (I)  **3)  /3. 

**4-X1  (IJ  **4J*  (Y  2  (IJ  -Y1  ( I),  J,/^ 

- r  :i ;  1 1  1  *  1 

T,°=P9S  (1)  *  { 
no=?P5  (I)  * 
p11=PPS  (3  * 

v  2  ]i  (  *  *  2-  X 1  (IJ**2  J  *  (Y  2  (I )  *  *3- Y1{I  j 
(X2(I)-X1  (I) )  *  (Y2  (I)  **4- Y  1  (I)**4)/4. 

(X2  « I)  **4-X1  (IJ  **4)*  (Y2  (I)  **^-Yl  (I)**2)/8. 

ji=(nh-i 

L 1=  (  M  2-  1 
K  1=  f  NI-  3-  1 

T*3 TT" 

)  *341 
)  *3  ♦  1 

I  4  1  .  4 _ 

-j—z-  -A  r-^-3-7  ‘7  l  x  1  \ 

A2= A**2 
B2=E**2 

A  F=  A  *  P 

- — 

33=E**3 

A  E  2=  A  *  P  *  * 

J  -  J  T  ' - 

2 

A2B=  A**2*B 
AE3=A*3**3 
A3«=  A**3*R 

P  (J1)  =P  (J1)  +  (F1-0.75+FU/A2-0.  25*F5/AB-0.  75«*F6/B2+ 0.  2  5*F7/A3  + 

*  0.  375«t;8/A2E*0.3  7  5»F9/AB2*0.  2  5»F10/B3-0.  1 2 5*  FI  1/A3B-0.  125* 

*P  (J1  ♦l{=F(J1  +1)  ♦  (-F30.  5*  F5/A  +  F6/B-  0.  5* F9/A B- 0. 2 5*F 1 0/B2+ 

*  0.  125*F  I  2/ AB2) 

P  |.lUil=PMU;i  +  (F2-5U/A-0.  5*F5/B+0. 25*F7/A2*C.  5*F3/AR- 

-* - 0.  4?5*Pl1/A2D) - 

P  ( L  1 )  ■=  P  (11)  +  (0. 25*F5/AD+  0.  75*F6/B2-0.  375*F9/A2B-0. 3 75*F9/AB2- 

*  0. 2S*F10/F2+0. 125*F1 1/A3B+0. 125*F12/AE3)  . 

PfLl-*  1)=P  Ll+1)  MO.  5*  F  6/5-0.  25»F9/A3-0.  25»?10/B2+0.  125»F12/AB2) 

P(KlI=P  JKUMO.  752jt/A2^o72’5*F5%B^of55*?7/A3-0^3^5*F8/A2B- 
0 .  37  5*>'9/AE2  +  0.  125*F1 1/A3B*0.  125*F  12/ AB3) 

*  '  «  c.T <r  .,.n  c.t.fty.n_n  .  1  51*?  . 

.  _  J  1 1/A2B) 
125* 


T  v/  f  -  T  f  7/  ut  <  t  w,  1 2.  j  T  :  i  i/<waa»u#  i  il/  noji 

P(Kl+1)=F (K  1  ♦  1 ) ♦ (-0. 5*F5/A+0. 5*F9/AE-0.  125*F12/AB2) 

— Pir-K?Ur]Ki*2{  rF-ft.  5»rV*>0. 25»r7/A2»0.  25»F8/AB-0.  1 25*  F 
P  i:i1)=P  (SI)  ♦  (~0.25*F5/AB+0.  37 5* F 0/A2B+-0.  375*F9/AB2~0-  125 
*  FI  1/A3P-0. 125*ri2/AB3) 

P  (Hi  ♦1)  =  F  ('ll  ♦  1)  ♦  (0. 25*F9/AB-n  ■•■>  ***  i  *»  ,» 


0.  1 2  5*F  1 


2006  CC  NT  I NUI 
EFTURF 
FND 


2/AB2) 

4  <  f  ft  *1  T1  1 

»  1/ A  luj 
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- DO-3-Wt*-J-  Tl — 1 - 

.1000  P  (I  ♦  1)  =F  ( I  +  L  ) - E-  (I)  *SU  (J  +  L) 
I?=I2  +  1 
114=1  1-  1 

- ne— ^e-ei-3-=T2-;  rea-s - 

o=  (i-i)  *icn7-n 
f  ( x)  =  r  (i)/su  (j) 

_ 13=11-1 

3*103  no  30on  i=i,  in' 

V  (1+  L)  =  P  (1  +  1)  -F  (I)  «SL’  (J  +  L) 
13  =  1 3-  1 

- I-f— f -*■  -I  ■?)— 9-3-e*7-3frfr*7-3  00  5 - 

3  004  CON  .  I K  U  r 
3C05  14=14-1 
3C01  CONTINUE 


I=IEHP-!U1 
J=  (1-1)  *1D*17  42 
CO  3006  K=1,  1 1 

po7 = K  | 
cm  i  ;:ur 
RETURN 


♦ST (J+K- 1) 


SIR'S*  *IPM.R«.N2X.R34ftfIsW5.SHl«STai|  IW2 

1  ST  FE2,P#TDJl8,X,y,V,U<  NLAY  F  R  ,  1  1  ,  E2\  T1, 12, COUP,  Tl, 

2  XS1AB.X  ,YF,SlX,DSX,riX.DJHX.nPHX.FyX-LTrX,AGSX, 

3  - ffilTtaH  ,  tl  YVi;d frYvCW-'J  r?  YY  ,  LTDYrA-QfrY  ,obii  ,-r^i-T) - 

IKTFG*F:  CCFP 

D13INSICN  ST  F  1  (1DI19 , 3)  ,  ST  I  2  {1EJ19 . 3)  , SI RT  1  (I  DM  9,5)  ,STRB1  (IEM9,S) 
2  STPT2  (IDM9  ,  5)  ,5TRB2  (IDK9, 5)  ,  P  fIDS6),X  (IEH9)  ,  Y  (IEM9)  , 

6 - Vf2)  ,'J  (2)  ,  f  ^4^19)  ,F2  (kD?l9H^1  (ID“9f  - — 

TT  (IEM9)  #XF(1D«9)  ,YF  (IEM9) 

N1XY=F1X*NY 
F  2XY=N2X*  NY* F  1XY 

bo— u-3-e-e— i=ivi-Ei;9 - - - 

IF  fI.lF.H1XY)  GO  TO  4001 
IF  (I.LF.N2XY)  GO  TO  4002 
I1=VY 


f  '  GC 

i  TO 

i  4 0C3 

! 

P  12 

"  13 

14 

- oe 

!=  N  2  X 
l=MXY 
=  KY-K2Y 
i  TO' 40( 

^3 - 

4001  I1=NY 
I2=N  IX 


4  003  D 
S 

4  004  S 


IF  (I.ifl.II* 1.CR.I. EQ.I1 *13.  CR.I.EQ.  II*  (12-  1)  +  1-*I3. OR. I. EQ. I  1*12* 
13)  C-C  TO  4006 

IF  (7  .  FQ.  I  3+i  4.CF.1.IC-I3  4I4+1.0R.I.EC.I1*  (12-1)  +N1Y  +  I3. 

- rk.-Ti  2-1)  +  HtY-»l3+  1)  GO  TC  40Q6  — — - - 

I.FC-  (I-I3- 1) /1 1*1 1*2  1*13)  GO  TC  4007 


IF  (I.fC.  (I-I3-1)/I1*I1*1  1*13)  GO  TC  4007 
If  (I.FC-  (I- 1  3- 1) /l  1  *11+ I  1*1  3-N2  Y)  GO  TC  4007 
IF  (1.GI..I1*  II  2-  1  1  *13)  GO  TO  4008 


•'  I"  ■*  / 

B=  (Y  (2  ♦  1)  —  Y  (  )  )/ 2. 

A  1=  A/F 

31=E/A 


J=I*  i 
L=I 1*3 

STH1  (I,1L=16«  (B1*V(1)  *A1)_*P  (J-2)-8*  A*V  (  1)  *P  (J- 1^* 8*B*P  ( J)  -6*  A 1* 


r.S|i:}USl}JiP5n-(S:5,,:!JilHi5n-»iSf?SN5lJU«fU«** 


r.-B2  |1, 3)  =  0(5)  *  (-2*P  (J-2)  *<4*8*?  (J-  1)-4*A*P  (J)  ♦  2*  P 
1  r(J*3)*2*r(J*L-2)-4*p*F(J*L-1)-2*P(J*L+1))/AB4 

4  009  IF  (I.EQ.  {1-1 J-  1)/I1*1  1+ 1*13)  GO  TO  4008 

IF  (1.50.  (I-l'?-1)/I1*im*13*NY-M2Y)  GO  TO  4008 
4007  A  =  (X  (I  ♦  !  1)  -X  (I)  )  /2. 

E-  (Y  (I )  -  Y  (I-  1)  )/2. 


1  P 

2  P 
STFI 


P  (.1*1 


4010  IF 


4  C  09  A  =  (X  (I)  -x  (I-I1 
B=  (Y  (14  1)-Y  (I) 
A 1=  A  /P 


At 4= 4*A*E 

0=1*3 

1=11*3 


)/Al4 


=$T B 1  (1,3)  *  U  (1)  *  J-2*P  (J-L-2)  *4*B*P  JJ-L-1)  *2*P (J-L*1 ) *2* 
—  4* P*P  (J-  1  }~4*A*P  (J)-2*P  (J  +  1)  +  4*A*P  (J*3)  J/AB4 


STR2  (1,1 )  =ST  h  2  (I ,  1)  +  (-(.*B1*P  (J-L-2) -4*B*P  (J-L)  *6*  (Bl+V  (2)  *A1) 

1  P  (J-2)~ti*A*V  (2)  *P  (J-1)-8*B*P  (J)-6*A1*V  (2)  *P(JO)  -4*A*V  (2)  * 

2  P  (J*  2)  J/AP4 


— ^**.  ****** 


U01<4  IF  < 
H  1=  r 


..  ,  .lUll-K  . FC. 2)  GO  TO  14  0  16 

H1=f  1  (I)  *T1  (I)  **3/(12.*  (1.-V  (1 )  ♦  *2)  ) 
DO  UD17  0=1.2  , 


—  n  i  ’ w. t ji n  i 

=- SI BT 1  (1,0) . 

=  C. 

=  0. 


IF  (CC IP. FC. 1)  GO  TO  U019 

H1=i  1  (I)  *T  l  (I)  **3/  (12.*  (1.-V  (1)**2)  ) 

H2=F2  (I)  *T 2(1)  **3/(12.*  (1.-V  (2)**2)  ) 


STRT1  (I,J)=fil*C  .  *ST3  1  (1,0)/  (C*T  1  (1)  **2) 
STSDI  (I»J)=-STFI1  (1,0) 

STBT2  (I,0)=h:*6.*STB2  (I,0)/(C*T2  (1)  **2) 


H2=F2(I)/(12.* 

DO  H021  0=1,3 

STRT1  (1,.1)=H1*6.*STR1  (1,0)*  (T1  (I )  +2.  *ENA  1)  /C 


T2  (1) -2.  *ERA2)/C 
-T2(I)-2.  *EKA2)  /C 

.OB.I.EQ.  (I-I3-1)/I1*I1*I1*I3)  GC  TO 


* 

0-9i- 

IF  (I, 

GO  T 0 

4073 

1*1 3+ VY-  N2Y.OR.1.EO.  (1- 3  3-1 ) /II  *3  1 +11  + 

- 

mm 

GC 

n 

SI 

!  TO 
RT  1  ( 

pri  i 

4025  . -  —  —  i-  ^  ™ 

I.  1)=0. 

II  .11=0. 

o  i  c  IX  I  *  /  —  \s 

3TRR2  (I,  1)  =0. 


“025  IF  1)  .SQ.O.  .CB.STRT1  (I.  .EQ.O.)  GO  TO  4026 


4026  3TET  1  (1 , 3)  =  0. 
RTPM  (1 . 31  =  0. 
SIR!  2  (1 , 3J  =0  . 
STEP  2  (I  ,  3)  =0  . 


STRT1  (1,4)  =  (S1PT1JI,  1)+STET1  jI,2))/2.+ 

*  SORT  (0. 25*  (STET1  (I,  1)-STRT1  (1,2)  )**2+STRT1  (I,  3)* *2) 

STPP1  (I,4)  =  fiTKEI  (I,  1)  +S7RB1  (I,2)J/2.  ♦  _ 

STRT2  (I  ,  4)=  (  S^t  RT2  (1^1  )+^T  PT2^(I  ^2))/2i  2t~ 

■^577^?STllT3TiSr^3^Kn??!ilrJir,-'**2^Tlw21I^>'*2^ 

*  SQM  (0.  25*  (STET  1  (I,  1)  -STRT1  (1 , 2)  )  *  *  2+ STKT 1  (1,3)  **2 ) 
STEB1  <I,5)=  (STFB1  (I,  1) +STRB1  (1.2)  )/2.- 

- StTHTfC- 25  +  (ST R 3  1  (I,  1 )  -STffB  1(1  ;  2)  )*  +? +5T10BTf  *2 )“ 

STRT2  (1.5)*  (FTRT2  (1,1)  *5T»iT2  (I,2))>2.- 

*  SORT  (0.  25*  (STTT2  (1,1  j-STRT  2  (1,2)  )**2+STRT2  (1,3)  ** 2) 
3TPE2jIJ?5)  *  (SIRf2JI,  1)  ♦§TrB2jlt2)l/'2:- 


TT  (I)  =T1  (I)  +12  (I) 

4  COO  CCMINUT 

IF  (NSIAB. EQ. 1)  GO  TO  4100 


+  1 

J11=  (K1X-1)  *  1»  Y  ♦  1 
J33=  K  1 X*  EY 
cr.x=crx*cix 


OCX*  CCIX 
D5X= 1 2. *  EEX/ (EJWX**3*  (1.+FYXJ) 

r2x=6.*rrx/(uwx**2*  (i.+fyxj  * 


rc  4101  K=J1  1, J3 
J 1=  K 
J2=K  +  FY 


